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PACS 87. 15. A - Biomolecules: structure and physical properties - Theory, modeling, and com- 
puter simulation 

PACS 71.20.Rv - Electron density of states and band structure of crystalline solids - Polymers 
and organic compounds 

PACS 79.20.Hx - Electron and ion emission by liquids and solids; impact phenomena - Electron 
impact: secondary emission 

Abstract. - Radiation damage is an unavoidable process when performing structural investi- 
gations of biological macromolecules with X-rays. In crystallography this process can be limited 
through damage distribution in a crystal, while for single molecular imaging it can be outrun by 
employing short intense pulses. Secondary electron generation is crucial during damage formation 
and we present a study of urea, as model for biomatcrial. From first principles we calculate the 
band structure and energy loss function, and subsequently the inelastic electron cross section in 
urea. Using Molecular Dynamics simulations, we quantify the damage and study the magnitude 
and spatial extent of the electron cloud coming from an incident electron, as well as the dependence 
with initial energy. 
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Introduction. — Recent advances in the development 
of X-ray Free Electron Lasers (XFEL) offer a tantaliz- 
ing ability to do photon science using short intense X-ray 
pulses. As a consequence, a number of theoretical models 
describing the physics of extreme X-ray-material interac- 
tion have been presented [1-6]. The potential to use the 
XFEL to do 3D single bioparticlc imaging [1] has enhanced 
the efforts to theoretically describe the dynamics of a sam- 
ple being exposed to an XFEL pulse. On a time scale 
longer than 5 fs, most of the ionizations in a sample ex- 
posed to an X-ray pulse will not be due to primary photo 
ionizations, but due to inelastic electron scattering. In 
the aftermath of a single photo ionization, the photo elec- 
tron as well as consecutive Auger electrons will interact 
with outer shell electrons in the surrounding atoms, lead- 
ing to an electron cascade (electron shower) . This process 
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is exploited when detecting photons in photo- multipliers, 
but leads to undesired ionization and structure damage in 
diffractive imaging. It is therefore important to know the 
inelastic scattering properties and understand the dynam- 
ics of the electron generation when describing the damage 
process. The inelastic cross section provides information 
regarding the deposition of energy from an energetic elec- 
tron into the system. For free atoms and ions the elec- 
tron cross sections are known, and there are a few simple 
empirical expression presented in the literature that give 
good estimation of the free atomic cross sections, see for 
example rcf. [7]. 

Structural determination of biological relevant samples, 
such as proteins, viruses or cells, are of large interest in 
many branches of life science. In the near future, facili- 
ties such as the European XFEL source in Hamburg and 
the LINAC Coherent Light Source (LCLS) at the Stanford 
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Linear Accelerator Center, as well as table-top XFEL [8] 
have promised to provide short hard X-ray pulses, at wave- 
lengths suitable for atomic or near-atomic structural de- 
termination of biological molecules. To describe the sec- 
ondary electron production in biomolecular samples, we 
present a theoretical study of urea, in which the cross sec- 
tions were derived from first principle calculations, using 
a methodology developed and tested earlier [9, 10]. 

Method. — We perform simulations of the electron 
cascade in crystalline urea (CO(NH 2 )2) using the spatial 
electron dynamics program, ehole, that is a part of the 
GROMACS [11] Molecular Dynamics software package. To 
do these simulations we proceed according to the following 
work flow: (i) calculate the band structure and energy 
loss function (ELF), (ii) calculate the elastic and inelastic 
electron cross sections using the energy loss function for 
the latter, and finally (iii) simulate the electron cascade 
from a single electron at a certain energy. 

Band structure and energy loss function from first prin- 
ciple calculations. The method used to calculate the in- 
elastic cross section is based on the ability to calculate 
the energy loss function from first principle calculations. 
Performing such calculations requires a sample with a rel- 
atively few atoms in the unit cell. We have chosen to use 
urea as model for a biological sample due to three major 
reasons. Urea has a well known crystalline structure [12]. 
Urea has an atomic composition of biological character, 
containing light elements such as carbon, nitrogen, oxygen 
and hydrogen, and having a density close to the density 
of cellular constituents. Furthermore, the urea crystal is 
among the simplest crystalline organic materials known, 
with 16 atoms per unit cell in the tetragonal space group 

The band structure of an idealized urea crystal is de- 
termined from first principles using structural data found 
in the literature [13, 14]. We use the ABINIT first prin- 
ciple computer code [15] to generate electronic structure 
results within the Local Density Approximation (LDA). 
In our calculations, we use a 12x12x12 reciprocal space 
grid and a plane wave cutoff value of 680 eV for standard 
LDA pseudo-potentials. Additionally, a scissors operator 
of 1.46 cV is used to correct the underestimation of the 
LDA band gap to optical data [16]. 

Electron cross sections. The generation of inelas- 
tic electron cross sections from energy loss functions is 
well known, and several methods are found in the litera- 
ture. The method we use here is the so called Tanuma- 
Powell-Penn model (TPP), based on Bethe [17] and Lind- 
hard's [18] work. We refer the reader to references [19-22] 
for a detailed information about the model. In earlier 
publications we have, in addition to TPP, employed an al- 
ternative model developed by Ashley et al. [23-25] . In the 
present work we chose to use only the TPP approach, since 
Ashley model was found to underestimate the number of 
generated electrons [26] . The approach of using calculated 
energy loss functions when determining the inelastic cross 



sections and secondary electron cascades there from, has 
been presented and evaluated in the case of two semicon- 
ducting materials in an earlier work [10]. To our knowl- 
edge, experimentally determined electron loss functions 
for urea are not available and we are in this work obliged 
to rely on an ELF calculated from first principles for the 
generation of the electron cross sections. Estimations of 
the inelastic electron cross sections have been presented 
earlier, both based on calculated and experimentally de- 
rived optical constants for biological material [27] , as well 
as direct fitting to inelastic experimental data [28-30] . Pi- 
oneering work in determining the energy loss function in 
biomatcrial from first principle calculations was initiated 
by Zaider and coworkers [31]. Tanuma et al. [28] present a 
study of the electron mean free paths between 50 and 2000 
eV for 14 organic compounds, comparing and evaluating 
two fittings to the so called TPP approach (employed in 
the present work). 

The elastic cross section is calculated with the Barbi- 
eri/van Hove Phase Shift package using a partial wave 
expansion technique [32,33]. The elastic cross section has 
limited impact on the final characteristics of the gener- 
ated cascades. Further details on these calculations are 
presented in [9,26,34]. 

Electron generation and dynamics. Once we have de- 
termined the electron cross sections, we simulate the elec- 
tron cascade as follows: (i) an initial electron with a cer- 
tain energy is generated. For each time step we calculate 
the probability for inelastic or elastic collisions from the 
cross sections. The probability for any collision event is 
compared to a random number which determines whether 
the electron is scattered or not. (ii) When the electron 
scatters inelastically the energy loss of the electron is 
possible through two channels: (a) if the energy loss is 
higher than the binding energy of the outer shell electron 
(E\ ost > -E-bind), an electron carrying a kinetic energy equal 
to the energy lost by the first electron minus the binding 
energy (E kin = £i ost - E hind ) is released. Now the sys- 
tem contains two electrons with the energy of the initial 
electron minus the binding energy of the second electron. 
These two electrons are now separately treated as the ini- 
tial electron was. The cascades evolve independently for 
each electron, (b) If the energy loss by the initial elec- 
tron is less than the binding energy (i?i os t < -E"bind)> the 
initial electron changes its path and loses the energy to 
the system, (iii) In the case of an elastic scattering event, 
the initial electron keeps its kinetic energy and changes 
its path, according to the elastic differential cross section. 
The generation of new electrons through impact ionization 
carries on until all electrons have a kinetic energy lower 
than the binding energy. Each electron trajectory is ex- 
plicitly followed, and the kinetic energy and momentum 
of each individual electron is known at every time step. A 
time step of one attosecond was used in the simulations. 

For computational simplicity we have assumed that 
all energy losses in the scattering events are absorbed 
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by the system. The model does not include hole scat- 
tering, electron-hole recombinations, nor clectron-phonon 
coupling. Neglecting the scattering of the holes accounts 
for a lowering of the total number of the ionizations events. 
In diamond this has been shown to underestimate the to- 
tal number of electrons generated nearly by a factor of 
two compared to experiments [26,35]. Previous results for 
water, however, have shown good agreement with radioly- 
sis data from experiments, albeit the measurements were 
performed on longer time scales [36]. 

Electron-hole recombination - where electrons and 
holes rccombine and deposit energy into the system as 
phonons, photons, or excitations - is in principle part of 
the scattering dynamics. However, the recombination co- 
efficient for urea was not found in the literature, and since 
the coefficients for other organic semiconducors span many 
orders of magnitude (see e.g. refs. [37-39]), estimates from 
analogous compounds will most likely be gravely mislead- 
ing. We have chosen to ignore recombinations in our cal- 
culations, knowing it could potentially overestimate the 
number of ionizations, as electrons that would recombinc 
cannot generate more electrons. However, the impact on 
the results will be limited as long as the carrier density is 
moderate, which is the case for all the electron cascades 
generated in this study. 

In an earlier experimental study we have fol- 
lowed the formation of secondary electron cascades in 
single-crystalline diamond caused by femtosecond X-ray 
pulses [35]. The experimental results showed good agree- 
ment with a theoretical calculation based on the TPP 
model [26], when excluding electron-hole recombination. 
Furthermore, as shown in the next section, the large band 
gap calculated for urea 4.74 eV is close to that of dia- 
mond 5.46 eV. Since the density of crystalline urea is lower 
than that of diamond (pdiamond ~ 3.5 g/cm 3 , p urea ~ 1.37 
g/cm 3 ), we expect the electron-hole recombination in urea 
to be lower than in diamond. 

Estimations of the vibrational spectra and the clectron- 
phonon coupling strength of the urea crystal require an ac- 
curate evaluation of the forces acting on the molecule upon 
intermolecular displacements. In its solid phase, urea has 
a planar configuration with vibrational modes ranging be- 
tween 2.5 meV to 439 meV at 10 K [40]. Hence, the en- 
ergy losses are expected to be somewhat larger than for 
diamond [26], but neglected here. 

Dissociative Electron Attachment (DEA), i.e. fragmen- 
tation following electron capture at subionization energies, 
has been demonstrated to play an important role in radi- 
ation damage of biomolecular samples [41,42]. No exper- 
imental DEA cross sections for urea are known, and al- 
though cross sections for several biomolecules can be found 
in literature, the majority are for the gas phase. Exper- 
imental DEA cross sections for tetrahydrofuran [43] - a 
rudimentary analogue to the deoxyribose rings in DNA - 
and acetone [44], both in the condensed phase, are nearly 
zero below 6 and 8 eV, respectively, and have peaks that 
are almost two orders of magnitude less than our calcu- 
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Fig. 1: Band structure of urea obtained by Local Density Ap- 
proximation. The energy bands show small dispersion and a 
direct band gap of 4.74 is observed at the gamma point (com- 
pared to the experimental value 6.18 eV [16]). A scissors op- 
erator of 1.46 eV is used to correct the LDA band gap to ex- 
periment and shifts the results to the right. The symmetry 
points in reciprocal space, listed along the x-axis are G=(000), 
X=(0§0), M=(§±0), A=(|||), R=(0§§) and Z=(00§). 

latcd elastic cross sections in the low energy regime. We 
have thus concluded that omitting DEA from our calcula- 
tions will have no or limited effects on the results. 

Results and Discussion. 

Band structure and energy loss function. The ob- 
tained LDA band structure is illustrated in fig. [l] by the 
dispersion of eigenstates along high symmetry directions 
for the urea crystal in reciprocal space. The dispersion is 
found to be in good agreement with results found in the 
literature. First principle calculations on urea have been 
presented earlier [46,47]. There are some minor disagree- 
ments between the band structure presented in fig. [T] and 
those presented in ref. [46,47], however due to the integra- 
tion over the energy loss function with respect to energy 
(see ref. [34] for details), the impact of small differences in 
the band structure on the calculated inelastic cross section 
is negligible. 

A direct band gap of 4.74 eV is found at the gamma 
point (6.18 eV experimentally [16]), other parameters cal- 
culated are the Fermi energy with respect to the bottom 
of the valence band (Ef = 12.81 eV) and the valence band 
width (E v — 9.71 eV). The energy loss function for zero 
momentum transfer is to our knowledge a feature absent in 
the literature, we obtain it from a linear response analysis 
within the random phase approximation and it provides a 
detailed structure for the energetics of the crystal within 
low-energy regime. 

The calculated energy loss function for urea, 
Im[— e(q = 0, E)^ 1 } is presented in fig. [2] In the 
absence of experimental data for urea, we provide a 
calculated ELF for ice and compare it to experimental 
data (inset in fig. pj. First principle calculations of 
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Fig. 2: Calculated energy loss function (ELF) of urea. The 
inset compares the calculated and experimental ELF [45] for 
ice. Due to the fewer number of van der Waals interactions in 
crystalline urea as compared to ice, we expect the calculations 
for urea to show better agreement with experiments than ice 
does. The first plasmon peak in urea is observed at 7.65 eV. 

molecular crystals, where hydrogen bonds and van der 
Waals interactions are abundant, has been much debated 
in the past years without a consensus being reached. The 
debate regards which parametrization of the exchange 
and correlation contributions to the Coulomb potential 
is the most suitable, see for example reference [48] for a 
discussion on cohesive energies. Our LDA band gap and 
the characteristically flat dispersion of eigenstates is a 
good indication that the optical transition probabilities 
are being well described. Also, the close agreement be- 
tween our calculations and experimental data for ice [45] , 
as seen in inset in fig. [2] gives our urea results further 
support, as does a comparison with polystyrene [30] 
calculated based on a fit to optical data [49] . For higher- 
energies (>50 eV) the cross sections are generated within 
the Free Electron Gas approximation. 

Inelastic electron cross section. The electron inelas- 
tic cross section for urea, together with those of water 
and diamond are presented in fig. [3] The water and dia- 
mond cross sections are presented as a comparison (taken 
from ref . [9] ) , both based on experimental energy loss func- 
tions. Being slightly denser than water, however with sim- 
ilar electronic properties, we expect the inelastic electron 
cross section to be higher than that for water (fig. [3|. 
Comparing the inelastic electron cross section for urea to 
the estimated cross sections for other organic compounds, 
such as nucleic acids [29,30], we can also conclude that our 
calculations agree with earlier estimations of similar com- 
pounds. As compared to inelastic cross sections for free 
atoms, the cross sections for diamond (calculated using 
the approach employed in this work [34]) is around 50% 
lower than that for free carbon atoms [50] (experimentally 
determined) . 

Ionization and electron dynamics. In contrast to ear- 
lier calculations of secondary electrons in water [9] and 
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Fig. 3: Inelastic cross-section for electron scattering in urea. 
Diamond and water are plotted as a comparison, taken from 
ref. [9] and [34] respectively. 

diamond [34], urea generates more secondary electrons 
from an initiating electron with the same energy. An elec- 
tron with a kinetic energy of 278 eV generates slightly 
above 12 secondary electrons in urea, whereas in water 
and diamond only 9 and 6 electrons are generated re- 
spectively. The number of electrons generated from a 
single impact electron for a certain compound is known 
to be related to the band gap of the material [26, 51] as 
-^electron ~ -Emit /(C-Egap) , where C is a compound specific 
constant. Fig. [4] illustrates how the average energy in the 
electron cloud drops down and equilibrates within a few 
femtoseconds. The figure shows the evolution of the en- 
ergy distribution from an originating electron of 500 eV. 
The electron cascade in urea reaches energies below the 
ionization threshold after around 5 fs, which is twice as 
fast as seen in water. The electron energy loss is related 
to the density of the crystal, and urea is therefore expected 
to equilibrate faster than water. The main Auger decay 
lines in urea are at around 250 eV (carbon) , 400 eV (nitro- 
gen) and 500 eV (oxygen). For these energies, the electron 
cascade generated by the Auger electron comprises of ap- 
proximately 12 electrons (carbon), 18 electrons (nitrogen) 
and 24 electrons (oxygen). 

In a small biomolecular sample exposed to an X-ray 
beam, the photoclcctron is likely to escape the system, 
and it is the Auger electron only that will generate the 
secondary electron cascade. In such systems the sample 
will generate approximately as many electrons as a pure 
water sample, as water has slightly lower average Auger 
energy and a slightly lower density. This makes water a 
good model system for simulating X-ray damage as was 
done by Bergh et al. [2]. For larger systems however, the 
photoelectron does not escape, and using water as a model 
system underestimates the number of generated electrons 
by about 30%. 

In the energy regime that is investigated in the present 
work, we found that the total number of electrons gener- 
ated after they have attenuated to subionization energies 
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Fig. 4: Time evolution of the energy distribution in secondary 
electron cascade, generated from an initiating electron with a 
kinetic energy of 500 eV. The total number of electrons is nor- 
malized to 1 for every timestep (by dividing by the total num- 
ber of electrons in the system). The figure shows an average 
over 100 simulations. 

is linearly dependent on the energy of the initiating elec- 
tron, as follows 

M — ^ init f 1 ^ 

Jv electron . r I- 1 - 1 

4.5-Eg a p 

Similar behaviour has been shown before for semiconduc- 
tors [10,51,52] as well as for diamond [26]. Damage due 
to ionization is related to the density of ionization events. 
The radius of gyration (i.e. the average distance of the 
electrons to the center of mass) of the generated electron 
cloud is not linear dependent on the initiating electron 
energy, but rather as a second order polynomial of the en- 
ergy (fig. [5]). This implies that the electron cloud from a 
low energy electron is much more dense, than that from 
a high energy electron, and therefore more damaging, de- 
spite the lower number of total electrons generated. The 
size of the ionizing electron cloud is quantified by the ra- 
dius of gyration of the cloud. This relates to the energy 
of the initiating electron, which serves as a basis for un- 
derstanding how the damage process depends with size in 
various crystals. 

Since the cross sections are calculated for a neutral sys- 
tem, we do not explicitly create positive charges in the 
system and only account for the number of electrons lost. 
For a highly ionized system, the use of atomic cross sec- 
tions will no longer be accurate, a plasma approximation 
is then needed [6] . This lapse in the model is likely to have 
an impact on the total number of electrons generated and 
remains to be studied. Furthermore, omitting the posi- 
tive charged ions will lead to a larger electron cloud than 
expected, since we cannot account for shielding effects. 
However, for a system that is only fragmentary ionized, 
the impact of disregarding the positive charges should be 
negligible. 
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Fig. 5: Radius of gyration (average distance from center of 
mass) of the electron cloud as a function of the energy of the 
initiating electron, after 50 fs. Each dot represents the average 
over 50 simulations, the line is a fit to the dots, a second order 
polynomial: y(x) = 0.18 + 0.013a; + 7.14 x 10" 6 :r 2 . 

Conclusion. — We have calculated the band structure 
and energy loss function for urea crystal from first princi- 
ples. Employing a methodology presented and tested ear- 
lier, we have determined the inelastic electron cross section 
and estimated the number of secondary electrons gener- 
ated from an incident electron at different energies. We 
are confident that the similarity of the ELF with experi- 
mental data and previous calculations on similar molecules 
implies our predictions are good models of reality. Our 
simulations show that urea generates 30% more electrons 
than water, from an initiating electron carrying the same 
energy. Such a difference is expected considering the den- 
sity difference of urea and water. Urea bears more simi- 
larities with the biological matter (such as cells, viruses) 
in both elemental composition and density, making it a 
more suitable object of theoretical study than water or 
diamond, especially when addressing large systems. Fur- 
thermore we describe the size of the electron cloud gen- 
erated from a single electron as a function of the electron 
energy, showing the spatial extent of the radiation damage 
due to secondary ionization. Secondary electron genera- 
tion is an important effect for the damage caused on a 
biomolccular sample by an intense X-ray beam, and it is 
therefore necessary to understand these processes in order 
to predict the damage and aid the structural determina- 
tion of single molecules or crystals. 
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